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Abstract 

The importance of dipole penguin operator, On = {9s/8TT 2 )mi > sa l j iU RT a bG^ 1 ' , 
on B decays is demonstrated. It is shown that branching ratios for the inclu- 
sive decay b — ► sss, the semi-inclusive decay b — > scf>, and the exclusive decay 
B — > 7TK are enhanced by 20% to 30%. A useful consequence of this calcu- 
lation is the inclusive b — * sss spectrum which will be essential in isolating 
pure penguin process. 
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Flavor changing one loop processes mediated by a gluon play an important role in un- 
derstanding the Standard Model (SM) [JJ-0]. It is by now well known that substantial 
corrections arise in these estimates from inclusion of the full electroweak effects (the 
7, Z penguins and the "box" contributions). Naively, one would think that in the B system 
electroweak corrections to gluon exchange may be negligible. This is, in fact, incorrect. The 
large mass of the top quark increases the relative contributions of the electroweak correc- 
tions. Our recent estimates for b — > sip process suggests a reduction in the rate of 20%~ 
30% f|||. In this paper we study another contribution, which has been mostly neglected so 
far, to nonleptonic B decays. This is the contribution from the Dipole Penguin Operators 
(DPO). Similar operators have been shown to give non- negligible contribution to e'/e |||| 
and hyperon decays 0. We consider three different types of decays in increasing model de- 
pendence, the inclusive b — > sqq and b — > sss, the semi-inclusive b — > s(f), and the exclusive 
B~ — > K°7c~ decyas. We find that the DPO gives important additional contributions to 
these decays and enhances their branching ratios by up to 30%. Any future quantitative 
analysis of B decays that relies on full electroweak Hamiltonian will have to include the 
DPO contributions. An important consequence of this calculation is the inclusive b — > sss 
spectrum which will be useful in isolating pure penguin process. 
The effective Hamiltonian 

The effective Hamiltonian up to one loop level in electroweak interaction for charmless 

B decays can be written as 

AC 12 
#ab=i = -^-[V uh V*{ Cl O x + c 2 2 ) - V tb V* CiO^ + H.C. , (1) 

V2 i=3 

where Oi are defined as 

Ox = qalpLfpffffLba , 2 = qj^Lff'fLb , 

03(5) = ti »Lb rfin L { R )q' > Qm = qai^Lbp ^ q'^^R)^ , 

q> q> 

3 3 
07(9) = -T/^Lb^e^^Rmq' , Q 8 (io) = 2?«7AE e 9'M L )9l > 

q' q> 

Oxi = ^m b qa^RT a bG^ , Q 12 = ^m b qa^RbF^ , 



where L(R) = (1 =F 7s)/2, / can be u or c, q can be d or s, q' is summed over u, d, s, 
and c. a and /3 are the color indices. T a is the SU(3) generator with the normalization 
Tr(T a T b ) = 5 ab /2. G£ v and F^ v are the gluon and photon field strength, respectively. O2, 
0\ are the tree level and QCD corrected operators. O^-q are the strong gluon induced 
operators. O 7 _i are the electroweak penguin operators due to 7 and Z exchange, and 
"box" diagrams at loop level. Ou^2 are the DPOs. 

The Wilson Coefficients (WC) q at a particular scale /1 are obtained by first calculating 
the WCs at my/ scale and then using the renormalization group equation to evolve them 
to \i. We have carryed out this analysis using the next-to-leading order QCD corrected 
WCs following Ref 0. Using a s (m z ) = 0.118, a em (m z ) = 1/128, m t = 176 GeV and 



|i~m(, = 5 GeV, we obtain 



= -0.3125 , c 2 = 1.1502 , c 3 = 0.0174 , c 4 = -0.0373 
= 0.0104 , c 6 = -0.0459 , c 7 = 1.398 x 10" 5 , 
c 8 = 3.919 x 10" 4 , c 9 = -0.0103 , c 10 = 1.987 x 10~ 3 . 



Cl 

c 5 



(2) 



The DPO coefficients at the two loop level have the following values ||11||: 



en = -0.299 , C12 = -0.634 . (3) 

The above Hamiltonian allows one to study the full set of charmless nonleptonic B 
meson decays. It is clear that only in b — > s transitions, the DPO may be important 
because of the enhancement factor I^V^/V^V^J ~ 57 in Eq.(l) compared with the tree 
contribution. The three processes mentioned before are all of this type. There are two types 
of DPO contributions. One from On and the other from Oyi- The contribution from 0\i 
is suppressed by a factor of a em /a s compared with the contributions from On. We will 
therefore neglect its contribution in the following analysis. 
The inclusive decay: b — > sqq 

Let us first consider b sss. This is a pure penguin induced process and does not involve 
uncertainties accosiated with hadronic matrix elements. It serves as a clear indication of 



the importance of DPO. It also might be feasible to isolate this process experimentally by 
selecting energetic kaons. In this process the gluon in the DPO has to break up into a pair 
of ss. In our calculation we will use the four-quark operator generated by On, 

H u = i^^V tb V*c u m b sa, u RT a bqYT a qk» , (4) 

where k = p b — p s in the gluon momentum, and q = s for b — > sss. 

Since there are two s quarks in b — > sss decay, we have to anti-symmetrize the amplitude. 
We have 

G 

A(b(p b ) -> s 1 ( Pl )s 2 (p 2 )s( Ps )) = --j=V tb V t * s [A 12 sa»Lbs2^Ls 
+ B^s^Lbs^Rs + C^Lbps^Rsp 

+ £>i 2 (p 6 +Pi) M s 1 i2r o 6s 2 y t r°s - (1 <-> 2)] , (5) 



where 



OL s C\\mi ,11 1, 
A 12 = 4c 3 + 4c 4 - 2c 9 - 2c 10 - " b - ( w ~ - - 

In l\ c s t 

a s c u m 2 b 1 

-D12 = 4c 5 - C7 - 



C12 = -2(4c 6 -2c 8 + 



2tt iV c s ' 
a s c n ml 1. 



2tt f 



a s c n m fe 1 

D12 = , (6) 

7T S 

where s = (p b — Pi) 2 -, t — (p b — p 2 ) 2 , and N c = 3 is the number of colors. The coefficients 
(A 2 i, B 21 ,C 2 i, D 21 ) are obtained by exchanging s and t in the above equations. The decay 
width is given by 

r = — — 1 „ / dsdt\A\ 2 , (7) 
(27r) 3 32mfi 1 1 ' K ' 

where the integration limits for s are integrated Am 2 and (m b — m s ) 2 , and the limits for t 



are (m 2 — s + 3m;? ± y (1 — 4m 2 /s)((m 2 1 — s — m 2 ) 2 — 4sm 2 ))/2. The integrand |A| 2 is given 
by 
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\A\ 2 = ^yfl^l 2 {[4iV c 2 (|A 12 | 2 + |A 21 | 2 ) + 8N c Re(A 12 A*J](s + t)(m 2 b -s-t) 

+ iV 2 [(4|I? 12 | 2 + \C 21 \ 2 )t(m 2 b -t) + (4|5 21 | 2 + |C 12 | 2 ) S (m 2 - s)} 

+ 2(iV 2 - l)[\D 12 \H{m 2 h -s) + \D 21 \ 2 s(m 2 b - t)](m 2 -s-t) 
N 2 - 1 

- 2^^ J Re(,D 12J D 2 * 1 ) S t(m 2 -s-t) 

- 4N c [Re(B 12 C* 21 )t(m 2 b - t) + ite^B^Mmg " s )] 
+ 2(N 2 - l)m b [Re{2A 12 D* 21 - C 12 D* 21 )s 

+ Re(2D 12 A* 21 - D 12 C* 21 )t}(m 2 b - s - t) } . (8) 

Here we have averaged the initial spin and color factors, and also taken into account the 
symetrization factor 2!. 

In our numerical evaluation, we will use \V ts \ ~ \V cb \ = 0.04, a s {m b ) = 0.25, m b = 5GeV. 
Since s quarks hadronize into kaons, it seems appropriate to use the constituent mass. 
We use m s = 0.5 GeV. Using these numbers, we calculate the relative branching ratio 
R = T(b -> sss)/T(b -> cev e ). Without the DPO effects, R = 0.018. When DPO effects 
are included, R = 0.024. We see that the DPO enhances the branching ratio by about 30% 
and should not be ignored. This result improves upon the result without QCD corrections 



discussed in Ref. [12|] 



It is essential to know the distribution of events with respect to invariant masses for 
selecting the desired events against background. To this end we present two useful spectra. 
In Figures 1 and 2, we plot the normalized event distributions with respect to diquark 
invariant mass x = (pi + p 2 ) 2 /(m b — m s ) 2 and quark anti-quark invariant mass y = (pi + 
Ps) 2 /(n^b — Tn s ) 2 . Although the effect of DPO on the overall branching ratio is significant, 
we see that its effect on the normalized distributions is very small. 

We now consider b — > sqq with qq summed over uu, dd and ss. This process may be easier 
to measure experimentally because this process can be isolated by selecting single kaons 
having high meneta. In this calculation we do not need to antisymmetrize the amplitudes 
for b — > suu and b — > sdd. We again use constituent masses with m u = ma = 0.3 GeV and 



5 



m s = 0.5 GeV. We find that the effect of DPO becomes weaker. It enhances the branching 
ratio by about 10% compared with the branching ratio without DPO contribution. The event 
distribution (shown in Fig 2) with respect to the qq invariant mass y = (p q +Pq) 2 / (m b — m s ) 2 
is similar to b — > sss. The DPO effect on the ditribution is negligible. The ratio R = T(b — > 
sqq)/Y{b — > cez/ e ) now also depends on the b — > u tree operators (Oi j2 ) which involves the 
weak phase 7. We find R is between 0.11 ~ 0.6 when varying the angle 7 between 0° ~ 180°. 
The semi-inclusive decay: b — > s(f) 

From experience we know that processes involving hadronic matrix element are difficult to 
handle because there is no satisfactory method to calculate these matrix elements at present. 
To have an idea how important the DPO effects are, we use factorization approximation to 
evaluate the DPO contribution to b — > s<f>. 

In the factorization approximation, the effective Hamiltonian in Eq.(4) gives vanishing 
matrix element without Fierz transformation because ss operator is in color octect states. 
Only the Fierz transformed Hamiltonian contributes. After a Fierz transformation, we 
obtain 



H eff , F = -^^V tb V*^-—cnm b {2m b s a lfM Lsps a, ^Lbp, 



Gf a s y V * N Z- 1 
V2 47TA; 2 tb ts 

- Am b s a Rsps a 'Lb pi + (p£ + tf)(s a ^Lsps a ' Rbp, + s a Rs p s a ' '^Rbp, 

2 N 

- is a o, v Rs p s<*'YRhp. + is a YLs p s a 'a, v Rb pi )]{5% - j^ffT^') . (9) 

The second term in the last bracket gives vanishing contribution in the factorization approx- 
imation. We have 

A{p _» s0) = _^^ 2 V tbV *^^ Climbig ^ [mb sYLb + 4(1 + ^rKsRb 

m 2 b + 2ml-m m s m b 

- m s (l - 2 ^ js^Rb-A 7 ^—p'rsLb), (10) 

m 2 ml 

where e M is the polarization vector of 0. In the above, we have used the parametrization 
< 0|s7^s|O >= igfa, with g\ = 0.0586 GeV 4 ; and < <p\sa^s\0 >= 5(e^-pfc v ) with the 
quark model relation 5 = —2m s g^jm\. We assumed that the s and s quarks inside each 



have similar momentum p$/2. In this approximation, k 2 = (ml — mh /2 + m 2 s )/2. Using 
the approximation A(b —>■ s<f>) ~ A(B — > X s <ft), for N c = 3 we find that without the DPO 
contribution the branching ratio BR(B —>■ X s <p) is 1.0 x 10 -4 . With the DPO contribution, 
BR(B — > X s <p) = 1.2 x 10~ 4 . The branching ratio is enhanced by about 20% when the 
DPO contribution is included. We also checked the DPO effect with N c = 2 favored by 
phenomenological fitting of the experimental data []T3| | . In this case the branching ratio is 



about two times larger, but the DPO contribution only enhances the braching ratio by about 
10%. 

The exclusive decay: B — > K% 

Here we carry out a calculation for B~ — > tt^K as an example. We, again, use the 
factorization approximation to estimate the amplitude. We have 

A ll( B- «-K°) = -^wS^Pun* 
x [2m b < K°\s'y IM Ld\0 >< 7r~|g7 M L5|5~ > 
- Am b < K°\sRd\0 >< n-\dLb\B- > 
+ (Pb+Ps)„(< K°\s lfl Ld\0 >< n-\dRb\B- > 
+ < K°\sRd\0 >< n-\d^Rb\B- > 



i 



< K°\sa Ufl Rd\0 >< n'\d^ u Rb\B- > 



+ i < K°\s-f u Ld\0 >< Ti-\da v ^Rb\B- >) . (11) 

In the above we have neglected the annihilation contributions of the form < 7r~.K' |,sr 1 'u|0 >< 
0|wr2|&|i?~ >, where Tj indicate the appropriate gamma matrices with corresponding 
Lorentz indices. These contributions are generally much smaller than the ones considered 
here M. 

We proceed with the Lorentz decomposition of the typical hadronic matrix elements: 

< K° \s^ 5 d\ >= iq^fx , with f K = 160 MeV , 

< 7T- \d^(l ± 7b)6| B- >= (p B + p^f + + (p B - p»)r , 

< 7T- da^{l + 75 )& B- >= -%2h{f B pl - p B p») . (12) 



We use the Heavy Quark Effective Theory (HQET) relation Jl4 h = (f + — f~)/Am b to 
relate h to f ± . 

In terms of the decay constant fx and form factors / , h, the B~ — > K tt~ decay 
amplitude is given by 

An{B~ - K\-) = ^£^V tb V*c n m b l^f K [((m| - ^)/+ + m 2 K f^ 

x ( mb 4- m b m 2 K ^m 2 B + 2m 2 K — m^\ 

I 2 (m b - m d )(m s + m d ) 8(m b - m d ) J 



m 



K 



4(m s + md) 



(2m| - im^)/^ + | (m| + 2m^ - rr£) 



+ r^BTr ("Ib + m K + m t ~ 2m|m^ - 2m 2 B m 2 K - 2m 2 K m 2 \ ] . (13) 



In the above we have used the same approximation for the kinemetics as in the case of b — ► 
s(j). We assume that the two quarks s and d inside K° share the meson momentum equally, 
and the b quark carries all the B meson momentum. In this case k 2 = (m 2 B — m 2 K /2 + m 2 )/2 

Using the form factors evaluated in Ref. [[0J, we obtain the branching ratio BR(B~ 
7i~K°) to be 0.83 x 10~ 5 without the DPO contribution. When the DPO contribution is 
included, BR(B~ — > tt~K°) = 1.0 x 10 -5 . Again, the DPO enhances the branching ratio by 
about 20%. The dependence on N c is weak in this case. 

From the above discussions we clearly see that the DPO has important effects on the 
branching ratio of B decays. These effects enhances the branching ratios by 20% to 30%. 
All analyses which involve numerical values of the strong penguin contribution are affected. 
These analyses need to be re-examined. However, since the DPO has the same transforma- 
tion properties under SU (3) flavor (and Isospin) symmetry as the ordinary strong penguin 
operators, many of the analyses related to CP violation in B decays which only depends on 
transformation properties under flavour symmetry |]6| JT6|Jl~7|] and processes where the strong 
penguin does not contribute B, will not be affected. 
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FIGURES 




FIG. 1. The solid and dashed lines are for the event distribution with respect to the normalized 
diquark invariant mass x with and without the dipole penguin operator contribution, respectively. 
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FIG. 2. The solid and dashed lines are for the event distribution with respect to normalized 
quark anti-quark invariant mass y for b — > sss with and without the dipole penguin operator 
contribution, respectively. The dotted line is for b — > sqq with dipole penguin operator contribution. 
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